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We introduce a dynamical model of a Bose-Einstein condensate based on the 2D Gross-Pitaevskii 
equation, in which the nonlinear coefficient is a function of radius. The model describes a situation 
with spatial modulation of the the negative atomic scattering length, via the Feshbach resonance 
controlled by a properly shaped magnetic of optical field. We focus on the configuration with 
the nonlinear coefficient different from zero in a circle or annulus, including the case of a narrow 
ring. Two-dimensional solitons are found in a numerical form, and also by means of a variational 
approximation; for an infinitely narrow ring, the soliton is found in an exact form. A stability 
region for the axisymmetric solitons is identified by means of numerical and analytical methods. In 
particular, if the nonlinearity is supported on the annulus, the upper stability border is determined 
by azimuthal perturbations; the stability region disappears if the ratio of the inner and outer radii 
of the annulus exceeds a critical value ~ 0.47. The model gives rise to bistability, as the stationary 
solitons coexist with stable breathers, whose stability region extends to higher values of the norm 
than that of the static solitons. The collapse threshold strongly increases with the radius of the 
inner hole of the annulus. Vortex solitons are found too, but they are unstable. 

PACS numbers: 03.75.Lm, 05.45.Yv, 42.65.Tg 



I. INTRODUCTION 



Matter-wave solitons have been created in Bose-Einstein condensates (BECs) in various effectively one-dimensional 
(ID) settings. First, these were dark solitons in repulsive condensates Then, bright solitons were created in an 
attractive BEC (lithium) 2]. This was followed by the making of gap solitons in a repulsive rubidium condensate 
loaded in a periodic potential, which was induced by the optical lattice (OL), i.e., interference pattern between two 
laser beams illuminating the medium 

A challenge to the experiment is creation of 2D matter- wave solitons. A natural problem in this case is the trend of 
solitons in multidimensional attractive condensates to be unstable because of the possibility of collapse in this setting 
Q. In theoretical works, several approaches were proposed to stabilize 2D solitons. One of them relies on the use 
of a full two-dimensional OL |^, or its low-dimensional (quasi- ID) counterpart 0, which can stabilize fundamental 
solitons. In addition, 2D lattices lend stability to vortical solitons including higher-order vortices, and "super- 
vortex" complexes 0; the latter are built as circular chains of compact vortices, with global vorticity imposed on 
top of the chain. Another theoretically elaborated approach relies upon the use of a nonlocal anisotropic nonlinearity 
induced by the long-range interactions between atoms with a magnetic momentum (chromium), polarized by an 
external field 

An alternative mechanism proposed for the stabilization of 2D matter-wave solitons is based on the Feshbach 
resonance (FR), which makes it possible to control the value of the scattering length, i.e., as a matter of fact, an 
effective nonlinear coefficient in the corresponding Gross-Pitaevskii equation (GPE), by means of an external magnetic 
field 0. Moreover, the FR may switch the sign of the nonlinearity (in particular, the FR-induced switch from repulsion 
to weak attraction was instrumental to the creation of bright solitons in lithium Q). Application of a low- frequency 
ac magnetic field may provide for periodic alternation of the nonlinearity sign in the GPE via the FR. It was predicted 
that the FR technique based on the ac field gives rise to novel states in the ID geometry and can stabilize 2D 
solitons, even in the absence of the external trap The same technique, if applied in combination with a quasi-lD 
OL potential, may also stabilize matter- wave solitons in the 3D geometry [T^ . 

It has been predicted and demonstrated in experiment that the FR can also be induced by a properly 
tuned optical field. Then, illuminating the condensate by two counterpropagating coherent laser beams, one can build 
an OL that will provide for periodic modulation of the nonlinearity coefficient along the respective spatial coordinate. 
Solitons in the corresponding one-dimensional GPE with the nonlinear OL were recently investigated in Ref. |15| . 
where stability regions for static solitons and breathers were found (motion of free solitons in the same model was 
recently studied in Ref. |0, and rigorous proofs concerning the stability of static solutions in this setting were 
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reported in work The soliton dynamics in the ID model with other configurations of the spatial modulation of 

the nonlinearity coefficient was studied in Refs. ^3 (unlike Ref. the nonlinearity coeSicient did not change its 
sign in the models considered in the latter works). 

Static spatial modulation of the nonlinearity through the FR, controlled by the properly shaped magnetic or optical 
field, may be tried as another means for the stabilization of 2D solitons, which is the subject of the present work. A 
natural form of am axisymmetric OL in the 2D geometry corresponds to the Bessel beam, i.e., a nondiffracting light 
signal in a bulk linear medium. In the case when the Bessel beam creates an effective linear potential in the equation 
of the GPE type with self-attraction, it has been shown that the corresponding radial lattice can readily stabilize 
various types of 2D solitons |0. However, our results show that, within a broad parameter region that we were able 
to explore, stabilization of 2D solitons by means of a nonlinear Bessel lattice, i.e., within the framework of the GPE 
whose nonlinear coefficient is g(r) = gQj„{ar), where r is the radial coordinate, go and a are constants, and J„ is the 
Bessel function with n = 0, 1, appears to be impossible - stationary axisymmetric soliton solutions can be easily 
constructed, but in simulations they all suffer either decay or collapse. 

Nevertheless, in this work we demonstrate that a simpler shape of the radial modulation of the nonlinearity, in 
which it takes a constant value, corresponding to self-attraction, inside a finite circle or annulus, and is zero (or 
corresponds to self-repulsion) outside this region, is able to stabilize axisymmetric 2D solitons. In addition to that, 
we will demonstrate that the model gives rise to bistability: the stationary solitons coexist with stable breathers, that 
feature persistent oscillations in the radial direction. In fact, the stability region of the breathers is larger than that 
of the static solitons, extending to higher values of the norm (number of atoms in the BEG). 

It should be said that, in the case of the nonlinearity controlled by the optical beam through the FR mechanism, 
the beam with the cross section in the form of a circle or annulus is not divergence-free, unlike its Bessel-shaped 
counterpart. However, this circumstance does not impede the physical realization of the model, as an effectively 
2D condensate can be easily trapped between two blue-detuned light sheets, which strongly repel the atoms, as 
demonstrated in the experiment |20|. The thickness of the corresponding "pancake" is a few microns, while its 
diameter is measured in hundreds of microns (at least) , hence the diffraction of the light beam within this range is 
completely negligible. 

The paper is organized as follows. In Section 2, we give the formulation of the model, and present numerical and 
analytical solutions for static solitons. The analytical part includes a variational approximation for the solutions in the 
general case, an exact solution for solitons supported by an infinitely narrow annulus carrying the nonlinearity, and 
predictions for the stability against radial perturbations, based on the Vakhitov-Kolokolov (VK) 2 1 . 22] criterion. The 
stability threshold for azimuthal perturbations is determined by a solution of the corresponding eigenvalue problem. 
An inference is that stability borders in the model with the nonlinearity supported on the circle are completely 
determined by radial perturbations, while in the annular model the upper stability border (in terms of the soliton's 
norm) is controlled by azimuthal perturbations. No stable solitons are possible if the annulus is relatively narrow, 
with the ratio of inner and outer radii exceeding a critical value ~ 0.47. In Section 3, we summarize results of direct 
numerical simulations of the stability of fundamental stationary solutions, which precisely confirm the existence of a 
well-defined stability region of the 2D solitons in the model's parameter space, predicted in Section 2. The bistability 
(coexistence of the stable stationary solitons and breathers) and the extended stability region for the breathers are 
also reported in Section 3. In Section 4, we briefly consider solitons with intrinsic vorticity, and conclude that all the 
vortices are unstable (the vortex splits in two fundamental solitons, each one then collapsing intrinsically). The paper 
is concluded by Section 5. 



II. STATIONARY SOLITONS 



A. The model and numerical solutions 



The GPE for the single-atom wave function ip in the normalized form is 

^^ = -^vV-5(0IV'PV', (1) 

with t time, the 2D Laplacian, and the nonlinearity coefficient shaped, by means of the external magnetic or 
optical field, as said above: 

9{r) = l}: P<r<R, (2) 

' \ 0, r < p 01 r > R. ^ ^ 
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FIG. 1: (a) Examples of stable soliton solutions with p = 0, ^ — —0.0399, and A*' = 5.59 (solid curve) and p — 0.5, 
p — —0.0648, and A'^ — 6.721 (dashed curve), (b) Chemical potential p vs. norm N for soliton families found numerically 
with p = 0, 0.2, and 0.5. For p — 0, the VK-stable portion of the solution, i.e., one with dN/dp < 0, is found in the interval 
N^^"'^ ^ 5.449 < iV < iVTowno. ~ 5.85. (c) p{N) curves predicted by the variational approximation for the same cases, p = 0, 
0.2, and 0.5. 

The number of atoms is determined by the norm of the wave function, 

noc 

N = 2Tr |V'(r)|Vdr. (3) 

Using the scaling invariance of Eq. we set R = 2, keeping p as a free parameter. Note that the model without 
the inner orifice, p = 0, is a universal one, as it contains no parameters. 

We also considered a model with the nonlinearity switched to self-repulsion, i.e., g(r) < 0, in the regions of r < p 
and r > R. However, we focus on the case with g — in these regions, as such a case is least favorable for the existence 
of solitons, hence it provides for results which are most relevant to the experimental realization of the scheme. 

Stationary solutions for fundamental solitons are looked for as ^ = 0(r)e~*'^*, with a real chemical potential /i, 
and a real function obeying the equation 

2p0 + (/."+ r^V + 25(r)03 ^0 (4) 

(the prime stands for d/dr). Equation Q is to be solved with the boundary conditions 0'(r = 0) = and 0(r = oo) = 
(the latter one implies that pL must be negative). The solution was searched for numerically by selecting the value of 
4>{r = 0) with which the boundary condition at r = cx) could be met. 

Two examples of the solution are displayed in Fig. 1(a), one for p = 0, i.e., the configuration with no inner "hole", 
and the other one with the "hole" corresponding to p = 0.5; in the latter case, the solution attains a maximum at 
r — p, having a shallow minimum at r = 0. Families of the soliton solutions are characterized by dependences p{N), 
which are displayed in Fig. 1(b) for p = and two nonzero values of p. These dependences predict a necessary 
stability condition as per the VK criterion dN/dpi < 0, i.e., parts of the solution families beneath the turning 

points in Fig. 1(b) may be stable (below, the turning point will be denoted as N = Nc°'"°'^^). In fact, the stability 
region exists due to the fact that the attractive nonlinearity acts in a finite region of space, r < R. 

In the absence of the inner orifice (p = 0), the soliton becomes narrow as p takes large negative values. In this case, 
the medium seems nearly uniform for the soliton, hence it approaches the shape of the well-known Townes soliton^ 
which is a universal weakly unstable localized solution of the 2D nonlinear Schrodinger (NLS) equation with the 
spatially uniform self-focusing nonlinearity [S^. Accordingly, the soliton's norm approaches the value A^xowncs ~ 5.85, 
which plays a critical role in the radial dynamics, being equal to the norm of the Townes soliton. 

B. Variational approximation 

The fundamental soliton solutions in the present model can also be obtained by means of the variational approxi- 
mation (see a review of the method in Ref. [23). To this end, we adopt the ansatz 



(5) 
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with an amplitude A and width w. The substitution of the ansatz in norm Q and Lagrangian of Eq. 



L = 2tt 



rdr, 



(6) 



yields iV = ttA'^w'^ (we use this relation to eliminate A in favor of N), and 



L = 2^iV - — 



1 



y2 27ru'^ 

Then, the variational equations, dL/dN — and dL/dw = 0, predict the following relations between the norm, width 
and chemical potential of the soliton. 



J 

N 



N 



The IJ.{N) dependence, predicted by Eq. ((TJ, is shown in Fig. 1(c) for several values of p. It is consistent with the 
numerical results displayed in Fig. 1(b), although the variational approximation predicts somewhat larger values of 
TV. 



C. The narrow-ring model 



The simple ansatz (jSl cannot predict the shape of the solution with the local minimum at r = 0, such as the one 
shown in Fig. 1(a) for p 7^ 0. The minimum becomes deeper as the nonlinearity-supporting annulus narrows, which 
corresponds to {R — p) / R ^ 0. As a limit form, one can take the GPE with the (5-functional nonlinearity support. 



(8) 



[the coefficient in front of the (5-function is scaled to be 1, cf. Eq. ©J. By final rescaling, one can again set i? = 2 in 
Eq. (jSJ, as was done above in Eq. (QJ, so as to cast Eq. (jSJ in a parameter- free form. It is relevant to mention that 
a BEG configuration in the form of a narrow ring was recently created in the experiment by means of an accordingly 
shaped magnetic trap j23 |. 

In the present case, the stationary wave function (j){r) obeys a linear equation, 



d^d) 1 d(t> 
dr'^ r dr 



0, 



(9) 



which must be solved separately for r < R and r > R. The inner and outer solutions, one with (j)'{r = 0) = and the 
other vanishing at r — + cx3, are to be linked by the conditions of the continuity of (/)(r) and jump of (f>'{r) at r = R, 
which follows from Eq. 



(/)'(r = R + 0)- (t>'{r = R-0) = -2\ 
Appropriate solutions to Eq. ||5J) are 



= R)f 



■^)^A 



Io(V^r) /loi^^R) , r<R, 
Ka lV^r) /Ko{V^R) , r > R, 



where Iq and Kq are the modified Bessel and Hankel functions, A is a constant, and the continuity of 0(r) at r 
is provided automatically. The substitution of expressions (tTTll in Eq. l(Tn|) yields 



A' 



Ki{z) . Ii{z) 



Ko (z) lo (z) 



The norm Q of the exact solution given by Eqs. Hll|l and (|12ll can also be calculated in an explicit form: 



N = ttJ-- 
2 



Ki{z) , /i(z) 



Ko (z) lo (z) 



KUz) lUz) 
Kl{z) 



(10) 

(11) 

= R 

(12) 
(13) 



Figures 2(a) and (b) display, respectively, an example of the solution, and the p{N) dependence plotted as per the 
exact expression (|13ll . 
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FIG. 3: The stability diagram for the soliton solutions. In the region between the two borders, the stationary solitons are 
stable - simultaneously according to the VK criterion, i.e., against radial perturbations (above the lower border), and against 
azimuthal modulations (below the upper border). 



D. Stability diagram for stationary solitons 



Figure 2(b) shows the existence of solutions with dN/dfi < in the model with the radial (5- function, which may be 
stable according to the VK criterion. However, it can only guarantee the stability against radial perturbations that 
do not break the axial symmetry of the solutions. On the other hand, it is well known that axisymmetric ring-shaped 
states may be easily subject to instability against azimuthal perturbations (see, e.g., Refs. p^). 

To study the stability against angular modulations in the general case [with g{x) taken as per Eq. (0)], including 
the 5-functional limit, as in Eq. JH)), we take a perturbed solution as 



(r) + S<j)+{r)e 



— iXt-\-irn9 



- d<p-{r)e 



iX* t — imO 



(14) 



where 9 is the angular variable, m is an integer perturbation index, x is a perturbation eigenfrequency, with * standing 
for the complex conjugation (x may be complex |25j|'). and 50±(r) are components of the respective eigenfunction. In 
particular, the instability threshold may correspond to x = 0, then the eigenfunction has (50+ = 5(j)- = d(j)o{r), and 
the substitution of expression (|I4|I in Eq. (jS)) and subsequent linearization lead to an equation for the zero mode, 



1 / Id m \ -.7 



54>Q = 0. 



(15) 



The instability threshold is achieved when real /i, found as an eigenvalue of Eq. IjlSI) . coincides with the actual 
value of the chemical potential of the unperturbed solution (/)(r) . This way, the threshold was identified for the lowest 
azimuthal perturbation mode, with m = 1 (in direct simulations presented in the next section, instability was observed 
solely against the azimuthal modulations with m = 1). 

The result of the analysis is summarized in Fig. 3, in the form of a stability diagram in the {p, N) parameter plane. 
The upper dotted border is the critical curve for the azimuthal instability with m = 1, found as described above, 
while the lower dashed curve is the existence and stability border for the soliton solutions, which is identified as a set 
of turning points of the /i(A^) curves in Fig. 1. Soliton solutions satisfying the VK criterion, dN/dfi < 0, exist above 
the lower border. Below the upper border, they are stable against the m = 1 azimuthal disturbances, i.e., the solitons 
are expected to be completely stable between the two curves. This expectation was verified by direct simulations, see 
the next section. 

Note that the border of the azimuthal instability in Fig. 3 is located, for p = 0, at a value of N which is identical to 
■^Townes ~ 5.85, i.e., in the case of p = (no inner orifice), the thresholds for the collapse in the radial direction, and 
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for the breakup of the axial symmetry in the azimuthal direction, are identical. The coincidence of the two thresholds 
for p = can be explained. Indeed, differentiation of Eq. Q in r shows that, for given (j>{r), the function (/)'(r) solves 
the following linear equation: 



0' = -2.9'(r) iq^ir)f 



(16) 



If g' = 0, Eq. 1)16(1 exactly coincides with Eq. (|15l) for m = 1, hence the function (/)'(r) may be identified as the 
corresponding zero mode. Of course, when 17 is a function of x defined by Eq. |j2J), which means g'{x) = S{r — R), 
the term on the right-hand side of Eq. 1)1611 does not allow (^'(x) to be the zero mode; nevertheless, in the limit of 
N A^Townos, the soliton shrinks to a size much smaller than R, hence becomes vanishingly small, along 

with the above-mentioned term. Thus, in the limit of A'^ = -^Townos, the function (j)'{r) provides for a solution to Eq. 
(|15|l with TO = 1, making N = ^Vtowiics the threshold of instability to the azimuthal perturbations with m = 1. 

A notable feature of the stability diagram in Fig. 3 is that the lower and upper stability borders meet and close 
down the stability region at p = p,nax ~ 0.95, which means that the nonlinearity-carrying annulus with the ratio of 
the inner and outer radii exceeding the critical value, Pmax/-R ~ 0.47, cannot support stable solitons. This conclusion 
implies that solitons cannot be stable either in model |(HJ with the radial (5-function. Indeed, detailed consideration 
of that model reveals the region of the azimuthal stability at A^ > 11.0 and /i > —0.0116, which entirely belongs to 
the upper branch of the p{N) curve in Fig. 2(b), with dN/dp > 0, i.e., the region is VK-unstable. 



III. DIRECT SIMULATIONS 



To check the predictions for the stability of the solitons, and examine the evolution of unstable ones, we have 
performed direct 2D simulations by dint of the split-step Fourier method, employing a basis composed of 512 x 512 
modes. The size of the integration domain was i x i = 60 x 60, with the center of the circle or annulus set at point 
{x,y) = (L/2,L/2), and the timestep At = 0.005. 

The simulations have confirmed the stability of the solitons in the region between the lower and upper borders in 
Fig. 3, and instability outside of this region. Figure 4(a) displays an example of the time evolution of \ip{x, L/2)\ (i.e., 
the profile of the cross section through the central point along the x axis) in a perturbed stable soliton, for p = 0. On 
the other hand, Figs. 4(b) and (c) demonstrate that (for the same case of p = 0) unstable solitons suffer collapse. 

However, unstable solitons [ones belonging to the upper, VK-unstable, part of the p{N) curve in Fig. 1(b), with 
dN/dp, > 0] whose norm is taken below a critical value, A'c"''''"'^^ w 5.99 (for p = 0), which is higher than the norm 
A'towiics ~ 5.85 of the Townes soliton in the two-dimensional NLS equation, neither collapse nor decay into radiation 
(in the NLS equation, a pulse with A^ < A'rowncs is bound to decay in the 2D uniform space). Instead, the unstable 
soliton rearranges itself into a stable breather. Figures 5(a) and (b) display an example of the evolution of breathers. 
In the simulations, the breathers remain stable indefinitely long, their oscillations getting more regular as N decreases. 
The amplitude of the oscillations, which we define as the root-mean square of the variation of the soliton's amplitude, 
A{t) = \u{x = y = L/2,t)\, decreases with A^, and it vanishes at another critical value, Nc°^^'^^ ~ 5.449. Up to the 
numerical accuracy, the latter one is precisely the smallest value of A^ at which the stationary solitons exist for p = 0, 

see Fig. 1(b). Thus, A^ = Nc™'^'''' is not only the point of the merger of the VK-stable and VK-unstable branches of 
the solutions, but also the one at which the breathers merge into the static solitons. 

We stress that the existence of the stable axisymmetric breathers up to A'c^"'''"^'^-' w 5.99 does not contradict the 




FIG. 4: (a) The evolution of \ip{x, L/2)\ (central cross section) in a stable soliton for p = 0, p = —0.189 and A'^ = 5.62. (b) 
An example of collapse of an unstable soliton, for p = 0, p = —0.0179 and N = 6.175. (c) The time dependence of the field 
amplitude, i.e., maximum value of |</)(a;, j/)|, for the same case as in (b). 
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(a) (b) 




" t A{0) 

FIG. 5: (a) The time evolution of the cross-section profile, |i/>(a;, L/2)], of a breather, for p = 0, fi = —0.0399, and A'^ = 5.59. 
The initial amplitude of the soliton is A{0) = 0.6. (b) Evolution of the amplitude of the breathing soliton for the same case, 
(c) The amplitude of intrinsic oscillations of the breather as a function of A{0). 



(a) (b) (c) 
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FIG. 6: Snapshots of contour maps of \ip{x,y)\ for an azimuthally unstable soliton, taken at t — 5 (a), t = 100 (b) and t = 122 
(c). In this case, p = 0.55, p. = -0.0782 and N = 7.055. 



fact that the symmetry-breaking azimuthal instability occurs, for p = 0, at iV > -/Vtowiics ~ 5.85, as explained above. 
Indeed, the latter pertains to the angular instability of the static solitons, but not breathers. 

A noteworthy consequence of these results is the bistability: in the entire interval of values of the norm, 
5.449 < N < A^Towncs ~ 5.85, where stable stationary solitons are found (for p — 0), they coexist with breathers. 
On the other hand, in the adjacent interval, 5.85 < N < Ncr^^'^'^'^ « 5.99, only stable breathers are possible (and no 
stable objects exist for N > 5.99). 

Stable breathers and bistability were found for p > as well. We note that stable breathers were also found in the 
model based on the one-dimensional GPE with a nonlinear OL [i.e., the nonlinearity coefRcient modulated in space 
as cos(A:a;)] In the latter model, bistability was observed too, as the breathers exist at the same values of the 

norm at which stable stationary solitons are found. 

As said above, all the solitons which are stable against collapse in the model with p = 0, are stable too against 
the azimuthal perturbations. Actual instability against the azimuthal mode (|I4|1 with m — 1 occurs at p > 0. To 
study the azimuthal instability in direct simulations, we used an initial condition in the form of a stationary soliton 
subjected to a weak angular deformation. Figure 6 displays a typical example of the development of the azimuthal 
instability for p = 0.55. As a result, the soliton does not split into fragments, which is a generic result of the azimuthal 
instability of vortex-ring solitons in uniform media 0|, but rather shifts from the central point, {x,y) = (30,30), to 
a position centered at {x,y) « (29,30). Because the norm of the soliton exceeds iVxownesj it then develops intrinsic 
collapse at the new position, where the hole does not essentially affect its dynamics. The shift of the soliton off the 
center and subsequent collapse were found to be a generic outcome of the development of the azimuthal instability. 
This feature can be easily explained by the fact obvious in Fig. 3: all the solitons which are subject to the azimuthal 
instability have N > A^Townos, hence they should collapse after being displaced away from the hole. 

As said above, direct simulations corroborate the stability of the stationary solitons in the region between the two 
borders in Fig. 3. We illustrate this conclusion in Fig. 7, which displays the time evolution of the field amplitude 
(maximum value of \ip{x,y)\) for p = 0.2 and three different values of the norm. The first soliton, with N = 5.702, 
belongs to the stability region in Fig. 3, and it is seen to be stable indeed. Two other solitons, with N — 5.994 and 
N = 6.01, are azimuthally unstable, which eventually leads to the collapse (after the spontaneous off-center shift, as 
shown in Fig. 6). Note that, as = 5.994 is close to the border of the azimuthal instability, the respective instability 
development time is large. 



8 




50 100 150 200 250 300 
t 



FIG. 7: The field amplitude vs. time, for weakly perturbed solitons with p — 0.2 and (a) — 5.702, fi = —0.0758, (b) 
iV = 5.994, /i = -0.164, and (c) = 6.01, p = -0.169. 




FIG. 8: (a) A typical example of profile (j){r) for the vortex soliton with p — 0, S = 1, and A'^ — 23.9, /j, — —1.33. (b) The p{N) 
dependence for the vortex-soliton family with p = and S — 1. 



IV. VORTEX SOLITONS 



Besides the fundamental solitons considered above, Eq. ^ also gives rise to vortex solitons, in the form of 
-0 = (f>s{r)e~'^'^'^~^^^^ , with integer vorticity S and real function (f>{r) satisfying the equation [cf. Eq. 

+ - S^r-^cf>s + 2g(r)0| + 2^^^s = 0. (17) 

In particular, in the model with the radial (5- function, see Eq. (jSJ, the vortex solution can be found in an exact form, 
cf. Eqs. mi and (O: 



0s(r) = A 



IsiV^r) /IsiV^R), r<R, 
KsiV^r) /KsiV^R) , r>R, 



2 V 2 



Is+iix)+Is-iix) Ks+i{x)+Ks-i{x) 
Is {x) Ks {x) 



The norm of this solution can also be calculated in an analytical form. 

An example of a vortex soliton, and the dependence /i(iV) for these solutions, are displayed in Figs. 8(a) and (b), 
for p = and 5=1. The figures show that a part of the solution family has dN/dji < 0, hence it is stable against 
radial perturbations, pursuant to the VK criterion. 

Comparing Fig. 8(b) to Fig. 1(b), one observes that the norm of the vortices is much larger than the norm of the 
fundamental solitons, which suggest that the vortex soliton may break up into a set of fundamental ones (as said above, 
this is a typical outcome of the development of azimuthal instability of vortex solitons in uniform media j2j|). Indeed, 
further analysis demonstrates that the vortex solitons with S = 1 are unstable against azimuthal disturbances with 
TO = 2 [cf. Eq. (|14|l ]. An example, displayed in Fig. 9 for p — Q, shows that the instability splits the vortex into a set 
of two zero- vorticity solitons, each then collapsing intrinsically, as its norm exceeds the critical value, iVrownos ~ 5.85. 
Before the collapse, the soliton pair rotates in the counter-clockwise direction. No example of a stable vortex soliton 
was found in the model. 



V. CONCLUSION 



The purpose of the work was to investigate the two-dimensional Gross-Pitaevskii equation in which the attractive 
nonlinearity is limited to a finite region in the form of a circle or annulus, including the case of a narrow ring. In Bose- 
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FIG. 9: Instability of vortex solitons is illustrated by a set of three snapshots of the contour map of \tp{x, y)\ for a vortex with 
S = 1, p = 0, AT = 23.9 and fi = -1.33, taken at t = 50 (a), 85 (b), and 90 (c). 

Einstein condensates trapped between a pair of blue-detuned light sheets, this configuration can be implemented 
through the Feshbach resonance by means of a properly configured magnetic or optical field which controls the 
scattering length of collisions between atoms. Using numerical and analytical methods, we have found a stability 
region for axisymmetric fundamental (zero-vorticity) solitons in the model, which is impossible in the case of the 
spatially uniform nonlinearity. It is noteworthy that the stability borders of the solitons in the model with the 
nonlinearity supported on the circle are completely determined by radial perturbations, while in the annular model 
the upper stability border is set by azimuthal modulations. The stability is limited to relatively broad annuli, with 
the ratio of the inner and outer radii smaller than a critical value, Pmax/^ ~ 0.47. Moreover, the model gives rise to 
bistability, as the stationary solitons coexist with stable axisymmetric breathers. The stability region of the breathers 
extends, in terms of their norm, to values exceeding the critical value corresponding to the Townes soliton. The 
collapse threshold strongly increases with the radius of the inner hole. Vortex solitons were constructed too, but they 
are unstable. Essentially the same results were obtained also for a model in which, outside of the circle or annulus, 
the nonlinearity is not zero but rather repulsive (that case is not considered in the paper, as the configuration with 
the zero nonlinearity is the most challenging one, as concerns the stability of solitons). The results reported in this 
work suggest a straightforward possibility to create stable two-dimensional matter-wave solitons in Bose-Einstein 
condensates. 
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